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In this paper we present a new formalism to implement the nuclear particle-vibration coupling 
(PVC) model. The key issue is the proper treatment of the continuum, that is allowed by the 
coordinate space representation. Our formalism, based on the use of zero-range interactions like the 
Skyrme forces, is microscopic and fully self-consistent. We apply it to the case of neutron single- 
particle states in *°Ca, ^°**Pb and ^"O. The first two cases are meant to illustrate the comparison 
with the usual (i.e., discrete) PVC model. However, we stress that the present approach allows to 
calculate properly the effect of PVC on resonant states. We compare our results with those from 
' experiments in which the particle transfer in the continuum region has been attempted. The latter 

, case, namely ■^'*0, is chosen as an example of a weakly-bound system. Such a nucleus, being double- 

■ magic and not displaying collective low-lying vibrational excitations, is characterized by quite pure 

neutron single-particle states around the Fermi surface. 



I. INTRODUCTION 

The accurate description of the single-particle (s.p.) strength in atomic nuclei is, to a large extent, an open issue (for 
a recent discussion see, e.g., Ref. ^Ij). Whereas in light nuclei either ab-initio or shell model calculations are feasible, 
^ [ in the case of medium-heavy nuclei we miss a fully microscopic theory that is able to account for the experimental 
(— I ' findings. Modern self-consistent models (either based on the mean- field Hamiltonians or on some implementation 
' of Density Functional Theory) do not reproduce, as a rule, the level density around the Fermi surface. The reader 
^ i can see, as a recent example, the results shown in Ref. [2]. Moreover, the fragmentation of the s.p. strength is by 
^ ■ definition outside the framework of those models. 

, In the past decades, much emphasis has been put on the impact on the s.p. properties provided by the coupling with 
various collective nuclear motions. The basic ideas leading to particle-vibration coupling (PVC) models in spherical 
OO nuclei, or particle- rotation coupling models in deformed systems, have been discussed in textbooks [3,]. These couplings 
provide dynamical content to the standard shell model, in keeping with the fact that the average potential becomes 
nonlocal in time or, in other words, frequency- or energy-dependent. We will call self-energy, in what follows, the 
dynamical part of the mean potential arising from vibrational coupling. This contribution will be added to the static 
CN [ Hartree-Fock (HF) potential. In this way, one may be able to describe the fragmentation and the related spectroscopic 
factors of the s.p. states, their density (which is proportional to the effective mass m* near the Fermi energy), the 
s.p. spreading widths, and the imaginary component of the optical potential. 

In, e.g., the review article [4| one can find a detailed discussion about the points mentioned in the previous paragraph, 
together with the relevant equations and the results of many calculations performed in the 80 's for the single-particle 
[ strength (mainly in ^"^Pb). These calculations are mostly not self-consistent and it is hard to extract from them 
■ ■ quantitative conclusions because of the various approximations involved. Certainly, they all agree qualitatively in 
pointing out that PVC plays a decisive role to bring the density of levels near the Fermi energy in better agreement 
with experiment or, in other words, the effective mass m* close to the empirical value m* « m. 

This enhancement of the effective mass around the Fermi energy, as compared to the HF value, is only one example 
of a phenomenon that can be explained by assuming that single-particle and vibrational degrees of freedom are 
not independent. Other examples, although not treated in the current work, are worth to be mentioned in this 
Introduction. Several works have identified the exchairee of vibrational quanta (phonons) between particles as one 
important mechanism responsible for nuclear pairing [S, Q. In the approach of 0, and references therein, one 
aims at explaining the properties of superfluid nuclei by taking into account both the pairing induced by the phonon 
exchange and the self-energy mentioned above (cf. also the discussion in Refs. [1, We also remark that important 
developments are under way ai ming at implementing ab-initio calculation schemes for open shell nuclei, based on 
self-consistent Green's functions ^11] or on the unitary correlation operator method [T5|). Along the same line, more 
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complicated processes can be explained by starting from elementary single-particle and vibrational degrees of freedom, 
and treating their coupling within the framework of an appropriate field theory: the spreading width of nuclear giant 
resonances, or the anharmonicity of two-phonon states (to mention only a few examples). The development of such 
a general many-body perturbation theory scheme could not avoid, so far, to resort to various approximations. In 
particular most of the calculations have employed simple, phenomenological coupling Hamiltonians. 

Recently, in order to calculate the s.p. st reng th, microscopic PVC calculations have become available, either based 
on the nonrelativistic Skyrme Hamiltonian ^13] or on Relativistic Mean Field (RMF) parameterizations fT^HH. The 
results seem to be satisfactory, in a qualitative or semi-quantitative sense, as they point to an increase of the effective 
mass around the Fermi energy. The results are clearly sensitive to the collectivity of the low-lying phonons produced 
by the self-consistent calculations. It is still unclear whether the results will eventually be improved by a re-fitting of 
the effective interactions or by the inclusion of higher-order processes. 

One of the limitations of all the PVC models that have been introduced so far, lies in the fact that they discretize the 
s.p. continuum (clearly, this means that the description of the vibrations themselves relies on the same approximation). 
Although in Ref. [16] a scheme to calculate the self-energy in coordinate-space representation had been proposed, 
there is at present no available result for the s.p. strength (let alone more complex physical observables) that 
avoids the continuum discretization. Consequently, the goal of the present work is to introduce for the first time a 
consistent description of PVC with a proper treatment of the continuum. In order to achieve this, the coordinate 
space representation is used. The current work is based on previous experience on how to treat the continuum within 
the linear response theory or Random Phase Approximation (RPA) (see Refs. [134111 )• 

The outline of the paper is as follows. In Sec. |lll we describe our formalism starting from the general formulation 
and stressing the implementation of proper continuum treatment. The main goal of the section is to display the 
equations that we have implemented and solved, and discuss the s.p. level density. In Sec. IIIIl the results for our 
three nuclei of choice are presented and discussed; whenever possibile, they are compared with experimental data. 
Finally, we summarize the paper and draw our conclusions in Sec. IIVI Some details of the calculations are shown in 
a few Appendices. 



II. FORMALISM 



A. Dyson equation in coodinate space representation. 

The particle-vibration coupling (PVC) Hamiltonian 0, in coordinate space can be written as 

Hpvc = jdr Sp{r)K{r)'y^^iJ;\ra)iJ){ra). 

(1) 

The density variation operator 5p{r) = p{r) — {p{r)) (where the brackets denote the ground-state expectation value) 
in second quantized form is given by 

5p{r) = J2 [SPnx{r)rl^ + Sp*r,^{r)r^x\ , (2) 

nX 

where f and f „a are the creation and annihilation operators, respectively, of a phonon n having multipolarity A 
and Sp„\ is the corresponding transition density, whereas k is the residual force. 

If the total Hamiltonian is H — Hp + Hpyg, where the term Hq describes uncoupled s.p. states and vibrations, the 
many-body perturbation theory [2^ [2^ can be applied. In particular, we assume that Hq includes the HF Hamiltonian 
for the nucleons and the independent boson Hamiltonian for the phonons (based on their RPA energies). We treat the 
term Hpvc as a perturbation using the interaction picture. We define Green's functions in space-time representation 
and we apply standard tools like the Wick's theorem to obtain the Dyson equation in terms of the unperturbed HF 
Green's function Gq and the perturbed Green's functions G: 

G{rat, r'a't') = Goirat, r'cr't') + ^ I dtidt2 / dridr2Ga{rat, riaiti)E{riaiti,r2a2t2)G{r2cr2t2,r'a't'). (3) 



The HF Green's function satisfies (w — ho)Go = 1, where ho is the s.p. HF Hamiltonian. The self-energy function is 
defined by 

T,{riaiti,r2cr2t2) = K{ri)G{riaiti,r2<T2t2)K{r2)iR{riti,r2t2), (4) 
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FIG. 1: The Feynman diagram for the self-energy function which corresponds to Eq. (O, in the approximation in which Go 
replaces G. 

where Rjr ^U r -fM) is the RPA response function (or phonon propagator) in the space-time representation, and is 
defined by ^ 

iR{rt, r't') = {^RPA\T{8p{rt)5p{r't')}\-^RPA), (5) 

where T denotes the time-ordered product and the formula stresses that the phonons are defined using the RPA 
vacuum \'^rpa) since this is exactly the phonon vacuum. We also note that the use of the Wick's theorem in the 
derivation of Eq. (jS]) implies the use of the causal representation of the Green's functions G and Gq, as well as of 
the RPA response function R. The connection between the causal representation with the retarded and advanced 
representations is outlined in the Appendices, where the causal functions will be denoted by G*^, i?*^. This label will 
be omitted in the main text, where we shall only make use of the causal functions. 
The Fourier transform of Eq. ^ is given by 



G{ra,r'a';uj) = Go(rcr, rV; w) + ^ J Jdridr2Go{ra,riai;uj)T,{riai,r2<72;u})G{r2a2,r'a']Uj), (6) 

(Ti(T2 

while the Fourier transform of the self-energy reads 

/oo 7 / 

— K{r)G{ra, rV; uj - uj')K{r')iR{rr' ■ a;') (7) 

due to the convolution theorem. 

A self-consistent solution of the Dyson equation involves the iteration of the two previous equations until convergence 
is reached. In practice, this is almost never done. In our work, since we explore for the first time the proper continuum 
coupling, we limit ourselves to the first iteration by replacing G with Gq in Eq. ([7]). 

We restrict our investigation to spherical systems in which static pairing correlations vanish. By taking profit of 
the spherical symmetry, one can use partial wave expansions and arrive at 

S(rcr,rV;a;) = 2^^^ — ^E;j(rr'; w)^^ — ; , 

Ijm 

(8) 

G{r(T,r'(j';uj) = 2^^^-^ — -Gij(rr' ;uj)^ — ; , 

Ijm 

(9) 
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FIG. 2: The Feynman diagrams associated with the perturbative expansion of the Dyson equation [Eq. ([3]) or Eq. ([6])]. 

where yijm{f'<^) = [^(^^) ® Xi/2{'^)] j,^- The latter equation holds evidently for Go as well. The Dyson equation can 
then be written as 

Gij{rr';uj) = Gajj{rr']iu) + JJdridr2 G'o,/j(rri;i:j)E;j(rir2; w)G;j(r2r'; w). (10) 
Similarly, the RPA response function and the residual force can be represented as 

Ri^rr';.) = ^^X^ R^rr' ; .)^^ , (11) 
K.{r) = K(r), (12) 
and consequently the self-energy can be calculated by 

^ij{rr-u:)= }^ __Go,z,„(rr ;w - w )-^ii?L(rr ; w ). (13) 

In our work, we start from the HF Green's function and RPA response function (together with the residual force 
k) and we obtain the self-energy from Eq. (|13p . Then, we also solve numerically the Dyson equation in the form (|10l) : 
for every energy of interest this equation can be cast in matrix form with respect to r and r' and solved as 

G= (l-GoS)-iGo. (14) 

In this way, the perturbed Green's function G contains the PVC perturbation up to infinite order, in keeping with 
the fact that it can be expressed by the Feynman diagrams represented in Fig. [2] 

B. Implementation of the proper treatment of the continuum 

1. Continuum HF Green's function and continuum RPA response function 

As already mentioned, our goal is an implementation of PVC that treats the continuum properly. In the case of 
atomic nuclei, in particular when local functionals like those based on the Skyrme interaction are used, considerable 
efforts have been made in this direction as far as the HF-RPA formalism is concerned. Indeed, the Green's function 
RPA has been formulated with Skyrme forces, with or without [l^ US] the continuum; the first self-consistent con- 
tinuum calculations have been presented in Ref. [2l| . In this context, proper treatment of the continuum means that 
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the Schrodinger equation including the HF mean field can be solved at any positive energy with the correct boundary 
conditions and, based on this, an exact representation of the HF Green's function Gq can be obtained. 
This unperturbed HF Green's function can be written as 

G^,,(rr'-E) = -A_iiy(r<; i?)«,,(r>; £;), (15) 

where uij{r; E) and vij{r; E) are, respectively, the regular and irregular solutions of the radial HF equation at energy 
E, r> (r<) are the larger (smaller) between r and r', and W{u^v) is the Wronskian given by 

( t j:^\9vij{r;E) duij{r;E)\ 



2m* (r) 

as 



is the (radial-dependent) HF effective mass which is defined as usual, in terms of the Skyrme force parameters, 
Ti^ If/ 1 \ / 1 M _ 1 



2™*(.) 2m ■ 4 ^^1^1+ 2-ij+^2^1+ 2^2j|p(r)-- U,(l + 2.0+^.(1 + 2x,)}p,(.). (17) 

In order to take properly into account the continuum effects also for the RPA phonons that lie above the threshold, 
the RPA response function appearing in the self-energy function will be calculated sclf-consistently, using the same 
Skyrme Hamiltonian used to compute the mean field. The details of the continuum RPA calculation have been given 
in previous papers [3 US]- We simply recall that two-body spin orbit and Goulomb terms, as well as spin-dependent 
terms, are dropped in the residual interaction. It is also necessary to convert the continuum RPA response function 
into the causal function, because normally the linear response theory (RPA) is formulated in terms of the retarded 
functions. This point is further discussed in Appendix B. 



2. Contour integration in the complex energy plane. 

Formally, the equations that appear in Subsec. Ill Al are defined in a model space which does not have an upper 
bound: in fact, the integrals over energy extend in principle from — oo to -|-oo, and single-nucleon as well as phonon 
energies have only a natural lower bound. However, the self-energy function does not converge if the upper limit on 
w' [Eqs. ([7]) and (1131) ] is extended to infinity, in keeping with the well-known ultraviolet divergence associated with 
the zero-range character of Skyrme forces. To avoid this, one must introduce a cutoff -Ecut- 

In order to make sure that only states below that cutoff contribute to the integrals in Eq. ([7]) and P^ . one can 
use the following procedure. By considering the expression ([T3|) for the self-energy function, one notices that the 
integral receives contribution from the poles of the causal HF Green's function and the causal RPA response function. 
The positions of these discrete and continuum (i.e., branch-cut) poles in the complex energy plane are schematically 
shown in Fig. |3l The blue dots and line represent the poles of the RPA response function, while the black crosses 
and line represent the poles of the HF Green's function. In order to pick up correctly the contribution of the poles 
below the cutoff, we must replace the integral in Eq. (|13p by an integral over an appropriate contour path. 
We have adopted the rectangular integration path displayed in Fig. [Sj which is similar to that employed in Ref. [l3l . 
It extends between — -Ecut to -Ecut on the real axis, and from to -rj' on the imaginary axis. 

It can also be shown that in this way one can reproduce the correct spectral representation of the self-energy 
function [cf. Eq. (jD3|) ] in the limit of a discrete system. 



3. Single-parttcle level denstty 
The level density associated with the HF single-particle levels can be defined by using the HF Green's function Gq 

as 

±1 

Po,ij{^)^ — / dr ImGo,;j(rr,w), (18) 
Jo 

where R is the upper limit of the integration. For bound states at negative energies, one is guaranteed that the 
result is stable with respect to increasing R. The sign ± guarantees that the level density is positive, i.e., the sign 
-|-(-) refers to particle (hole) states. This is equivalent to the definition pojj{uj) = S{uj — ej^j ■) and is normalized 
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FIG. 3: (Color online) Contour path C for the integration on uj' in Eq. and Eq. ([TSj. The blue dots and lines represent 
the poles of the RPA response function. The black crosses and line represent the poles of the HE Green's function. For the 
parameters r), r\ and _Ecut see the text. 



to 1 for bound states. In the absence of any potential, pQ^ij(uj) reduces to the free particle level density Pfree.iji'-^), 

2 

obtained by replacing Go.ij with the Green's function for the free particle which satisfies {oj — 2^)G/ree = 1- Gfree 
can be calculated either numerically or analytically using the same definition of Eq. (1151) . but with the use of the 

wave function for the free particle. This free particle level density is pfree,iji^) oc ^ 27ryzj large value of uj. 
For R — > oo, pfree,ij divcrgcs. On the other hand, po,/j(<^) tends to pfree,ij for large values of uj. It is then useful to 
introduce a new level density pojj(uj) by subtracting pfree,iji'^) [25l |. that is, 



±1 



TT 



R 



Po,;i(w) = — / dr Im (Go, /j(rr, w) - Gpree.ijirr^uj).) (19) 



In this way, the dependence on R is eliminated also for positive energies (cf. Fig. [S]). For uj < 0, there is no 
contribution associated with the free particles. As mentioned above, pojj coincides with the usual definition of the 
level density for the single-particle levels [25| 

S{uj — e|j°^-) for bound states (w < 0), 

Po.iji^) = { ^"^^(0) (20) 
for positive energy states (uj > 0). 

TT duj 

In an analogous way, we define the perturbed (HF+PVC) level density by using the solution of the Dyson equation 
Gij as 

±1 f 

pij{uj) = — dr lm{Gij{rr,uj) - GFree,ij{rr,uj)) (21) 

The peaks of the perturbed level density provide renormalized single-particle energies which include the effect of the 
particle-phonon coupling. In fact, if this coupling is small one can expect a simple shift of the HF peaks. Otherwise, 
the s.p. strength can be quite fragmented: the associated widths reflect the basic decay mechanisms that are the 
nucleon decay (providing the so-called escape width, or FI) and the spreading into the complicated configurations 
made up with nucleons and vibrations (providing the spreading width T^). 



III. RESULTS 

We shall present results for three nuclei: ^°Ca, ^''^Pb and ^''O. The effective Skyrme interaction SLy5 [IB] is used 
to calculate the HF mean field. The calculation of the RPA response is carried out exactly as described in Ref. fl^l in 
the limit of no pairing, using all the terms of the SLy5 interaction, except for the two-body spin- dependent terms, the 
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spin-orbit terms and the Coulomb term in the residual p-h force. In the calculations, the angular momentum cutoff 
for the unoccupied continuum states is set at Icut — for '"^Ca, and Icut — 12?i for ^^O and ^°^Pb, respectively. The 
radial mesh size is Ar = 0.2 fm. The values of the parameters used in the contour integrations (see the discussion in 
lnB2)) are 77 = 0.2 MeV [cf. Eq. ^] and rj' = 2t] [cf. Eq. ([rU))]. 
There are a few important issues that we wish to stress: 

1. Due specifically to the zero-range character of the Skyrme interaction, the self-energy diverges logarithmically 
as a function of the maximum energy of the phonons, as it has mentioned above. The first steps towards a 
systematic renormalization procedure have only recently been started to be worked out 27^. In this work, we 
shall take the usual view that the important couplings are those associated with the collective low-lying states 
and giant resonances. We shall then include phonons associated with the multipolarities 2+, 3^, 4+ and 5^, and 
set an upper cutoff on the phonon energies given by -Ecut = 60 MeV because no strong peaks are present above 
this value in the calculated RPA strengths. The way in which this cutoff is implemented has been described in 
detail in III B 21 

2. In the present scheme, the price to be paid for the exact continuum treatment, is that one cannot discriminate 
between the inclusion of collective and non-collective phonons. Actually, the diagrams shown in Figs. [T] and [5] 
contain terms that violate the Pauli principle, and these terms are larger when the phonons are non-collective. 
In other words, one could expect that in an exact calculation the correction of the Pauli principle violation 
cancels, to a large extent, the contributions from non-collective phonons. Although this point has never been 
clarified in the available literature, to our knowledge, on a quantitative basis, in most of the cases the usual 
view has been to take into account only the coupling to collective states. In the work that we quoted already 
(the most similar to the present one), namely in the recent calculation of particle- vibration coupling in ^''Ca 
and ^°^Pb of Ref. [ij], only phonons exhausting at least 5% of the isoscalar or isovector non-energy- weighted 
sum rules have been taken into account. Therefore, we can expect that the effect of particle-vibration coupling 
is larger when we calculate it with the present method, as compared with Ref. pj)] . We shall come back to this 
issue below. 

3. The momentum-dependent part of the particle-hole interaction had previously been neglected in the calculation 
of the particle-vibration coupling; then, in Ref. :13] it was shown that its effect is important (at least in the 
case of the SLy5 interaction), and that it can be reasonably accounted for within the Landau-Migdal (LM) 
approximation, by choosing the Fermi momentum kp as 1.33 fm~^ (that is, at the value associated with the 
nuclear matter saturation density). The LM approximation will be adopted in the following, with the same 
value of kp- 

A. Results for '"'Ca 

The first essential steps of our work consist in the calculation of the HF spectrum and of the RPA strength functions. 
The results are illustrated respectively in Table U (HE single-particle spectrum) and in Fig. 2] (isoscalar and isovector 
RPA strength functions associated with the multipolarities 2+,3~,4+ and 5~). In Table HIl we give the theoretical 
energies and transition strengths for the low-lying collective 3~ and 5^ states, comparing them with available data. 



Nucleus hole states [MeV] particle states [MeV] 





-48.3 


1/1 


-9.7 


Ipf 


-35.0 


2p| 


-5.3 




-31.0 


2pi 


-3.1 


Idf 


-22.1 


1/f 


-1.3 


24 


-17.3 






Idf 


-15.2 







TABLE I: Neutron single-particle energies in ^"Ca. 



In Fig. Owe compare the level densities pa^ij, PFree,ij and po,;j defined above [cf. Eqs. (|18ll2ip ]. in case of neutrons 
and for the quantum number gg/2- The results are shown for different values of the upper limit of integration, namely 
R = 15, 20 and 25 fm. The main goal of the figure is to illustrate the effect of the removal of the free particle level 
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FIG. 4: (Color online) Isoscalar (IS) and isovector (IV) RPA strength functions in ^"Ca for the multipolarities 2''',3 ,4"*" and 
5^. The IS 2^ and IS 3~ strength functions are reduced by a factor 10. 



Nucleus J'' Energy B{Q^j^°) Energy B{Q^j^°) 



40Ca 3- 
5- 



Theory (RPA) 
Energy 
[MeV] [e^fm^-^] 

4.13 1.13 X 10* 



5.25 2.22 X 10' 



,6 



Experiment 

KQT 

[MeV] [e^fm^-^] 
3.74 1.84 X W 
4.49 



TABLE II: The theoretical values for the energy and the transition strength B of the low-lying isoscalar 3 and 5 states in 
''"Ca are compared with the experimental data, that are taken in the case of the 3~ state from Ref. [2^ and in the case of the 
5" state from Ref. [H. 



density, that has been formally introduced above (cf. IIIB 3[) . In fact, it can be seen in panels (a) and (b) that for each 
value of R the smooth tails of pFree,ij{^) and po.ij converge, for uj between 5 and 10 MeV, towards the asymptotic 

value ^ 27iV" iiidicated by the dashed lines and by the symbol f{R,uj) in the panels. In panel (c) we show the 
level density po^ij{u!) defined by Eq. ((T^ . The free level density is eliminated (and the dependence on R with it) so 
that one can clearly identify the pure gg/2 resonance in the continuum. The inclusion of PVC leads to a quite big 
fragmentation of the strength, as it can be seen in panels (d) and (e) of Fig. [51 panel (d) is meant to mainly show 
that the 3~ states are the most important to produce that fragmentation, whereas in panel (e) we illustrate the effect 
of the subtraction procedure on the perturbed level density. 

Since the parameter rj introduced in our definition of Green's functions and response functions is one of the numerical 
inputs of our calculations, we have carefully checked whether the results are sensitive to the choice of its value. In 
Fig. [HI we compare the densities pij obtained with our standard choice of the the smearing parameter 77 = 0.2 MeV, 
and with 77 = 0.1 MeV, again corresponding to *"Ca and the Skyrme set SLy5, in the case of the quantum numbers 
gg/2 and pi/2. It is quite reassuring that the structure displayed by the peaks of the level density does not depend 
on the chosen value of 77. In principle, we may expect that a dependence of this kind shows up in the discrete part of 
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FIG. 5: (Color online) Neutron level density in the case of the g9/2 quantum numbers for Ca. These quantum numbers 
are associated with a single-particle resonance in the HF mean field. In panels (a), (b) and (c) we show respectively the 
quantities po.ij, ppree.ij, and po.ij defined in Eqs. (|18l2ip . for three different values of the upper integration limit R, namely 
15 fm (red curve), 20 fm (green curve) and 25 fm (blue curve). The asymptotic expression for the free particle level density, 
f{R,Lj) = ^ 2-nViU shown in (a) and (b). Note that in (c) the curves corresponding to different values of R practically 
coincide, because the free particle contribution is removed from po,ij- In (d) we show the perturbed densities pij, obtained by 
solving the Dyson equation with the inclusion of the coupling either with all multipolarities (blue curve) or with only the 3~ 
(red curve). They are compared with the HF density po.ij (green curve, R=20 fm) already shown in (a). In (e) we show the 
corresponding perturbed densities pij, obtained by subtracting the free-particle level density. They are compared with the HF 
density po,!j already shown in (c). 

the spectrum and vanishes when the continuum coupling becomes dominant: this effect can be to some extent seen 
in the high-energy part (above « 5 MeV) of the upper panel of the figure. 

In Fig. [71 we show results for the single-particle level density pij in ^°Ca, associated with various quantum numbers. 
The unperturbed level density is shown by means of the black curve, and displays sharp peaks of equal heights at the 
HF energies. We compare in the figure the results obtained by taking into account the coupling with all multipolarities 
(blue curve), or with 3~ phonons only (red curve). The first qualitative remark is that for the states lying close to the 
Fermi energy, both in the case of hole states (2si/2, ^^5/2 s-nd ld3/2) and bound particle states (2p3/2, '2pi/2^ ^^7/2 
and lf5/2), the strength remains concentrated in a single peak, eventually acquiring a spectroscopic factor, and the 
quasiparticle picture maintains its validity. This is not true when we consider states either more far from the Fermi 
surface or in the continuum, that is, at energies where we expect the single-particle self-energy to become larger. 

The hole states are in the left part of Fig. [71 For the aforementioned 2s and Id states there is only a shift of 
the HF peak. Instead, in the case of lpi/2 and IP3/2 states the strength is damped over a broad interval. It is 
interesting to trace the origin of this fragmentation, by restricting the summation over the phonon multipolarity and 
the angular momentum of the intermediate single-particle states in Eq. (jl3p . and analyzing the contribution of specific 
configurations. 

To help the following discussion, we depict in Fig. [8l selected contributions to the particle [(a) and (b)] and hole 
[(c) and (d)] self-energy. Since our equations and numerical codes are written in the coordinate space, these different 
contributions cannot, strictly speaking, be singled out. However, at definite energies it may happen that only one is 
dominant. We use Fig. [51 to recall that the fragmentation of a hole state can only be caused by on-shell contributions 
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FIG. 6: (Color online) Dependence of the HF+PVC level density pij [cf. Eq. (|2ip ] on the smearing parameter rj, for the 
states (upper panel) and for the pi/2 states (lower panel). Only the coupling with 3~ phonons is taken the account. The red 
curves show the results obtained with the standard value rj = 0.2 MeV, while the blue curves have been obtained with ?7 = 0.1 
MeV. We also show the HF results (dashed curves, obtained with rj = 0.2 MeV). 



associated with the coupling with other hole-phonon configurations [panel (c)], while the coupling with particle states 
[panel (d)] can only produce an energy shift. 

The unperturbed (si/2)^^ strength shows two peaks, associated with the lsi/2 (eisi/2 ~ —50 MeV) and 2si/2 
(e2si/2 ~ ~20 MeV) single-particle states. Due to parity and angular momentum conservation, the only intermediate 
hole-phonon configurations S1/2 holes can couple to are (d3/2)~^ ^ 2+ and (d5/2)~^ ® 2+. These coupling can lead to 
a strong fragmentation of the (lsi/2)^^ strength. In fact, the energy differences 61^3/2 — eisi/2 = —15.2 -I- 48.3 = 33.1 
MeV and 61^5/2 — eisi/2 = —22.1 + 48.3 = 26.2 MeV (cf. TablelJ) are close to the centroid of the isovector quadrupole 
strength (i?ivGDR = 29.2 MeV, cf. Fig. IJ]). On the other hand, although the (1^3/2)"^ <8)2+ intermediate configuration 
could in principle contribute to the fragmentation of the (2si/2)~^ state, in this case the energy difference is low, namely 
eid3/2 — eisi/2 ~ —15.2 + 17.3 = 2.1 MeV. In **°Ca there is no low-lying quadrupole strength, and therefore the HF 
strength is shifted but not fragmented. In a similar way, one concludes that the d3/2 and d5/2 hole strength, that lies 
close to the Fermi energy, cannot be fragmented. The deeply bound pi/2 and P3/2 states can instead couple efficiently, 
respectively to the (ld5/2)~^ 3^ and to the (ld3/2)^^ ® 3~ configurations. The relevant energy differences lie in the 
range 9 ~ 20 MeV, where one finds substantial 3~ strength (cf. Fig. S]). The case of pi/2 is analyzed in more detail 
in Fig. [9l including only 3^ phonons. It is seen that the full calculation (top panel, left) is very similar to the result 
obtained including only (ld5/2)^^ ® 3^ configurations (top panel, right). Coupling only to the lowest phonons of each 
multipolarity (bottom panel, left) leads only to a modest energy shift, again caused essentially by the (ld5/2)~^ ® 3^ 
configuration (bottom panel, right). The remnant of this configuration is visible in the figures, at about -25 MeV. 

The results obtained for the particle states are depicted in the right part of Fig. [T] As already mentioned, those 
lying close to the Fermi energy (2p3/2, 2pi/2, lf7/2 and lf5/2) can be well described within the quasiparticle picture: 
the associated single-particle strength shows a well defined peak, which is shifted from the unperturbed (HF) position. 
The situation is quite different for the unbound particle states lg7/2 and lgg/2 which are strongly fragmented. For 
these states, we expect that our proper treatment of the continuum should be particularly important. The case of 
the lg9/2 orbital, which in the HF calculation is associated with a low-lying resonance lying at about 2.5 MeV, is 
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FIG. 7: (Color online) The single particle level density fiij for neutron states in '*°Ca defined by Eq. (|19p for HF and Eq. (|2ip 
for HF+PVC. The left panel refers to hole states, the right panel to particle states. The black curve represents the HF level 
density, while the red curve and the blue curve show level densities resulting from the coupling to phonons. In the case of the 
red curve, only 3" RPA phonons are taken into account, while in the case of the blue curve 2"^,3~,4''" and 5~ phonons are 
taken into account. 



analyzed in more detail in Fig. 1101 We consider only the coupling with 3~ phonons, since they produce most of the 
fragmentation. By comparing the left and the right top panels of Fig. 1101 one concludes that the strong fragmentation 
of the resonant level is caused by the coupling with several intermediate configurations [namely (p3/2)®3~, (f5/2)'X'3~ 
and (f7/2)®3~]. The two satellite peaks found at w —6 MeV and at w —2 MeV are produced by specific configurations, 
associated with the lowest 3J~ collective state (cf. bottom panel, right). This phonon is responsible for about half of 
the total width (compare bottom panel, left with top panel, left). 

In Fig. [TT] we compare the position of the seven HF energy levels lying close to the Fermi energy with the position 
of the shifted levels, deduced from Fig. [T] Only for these levels a centroid energy is quite meaningful because, as we 
have already emphasized, for these levels essentially only one peak exists when one looks at the PVC results and the 
quasiparticle picture holds. Also for these levels, and for them only, one can attempt a comparison with the results 
of Ref. [Tsl that have been obtained through second-order perturbation theory and in a very similar scheme. The 
results are in overall agreement with those of Ref. [l^, although the magnitude of the present energy shifts is larger. 
In fact, while in Ref. [Ij] the shifts are typically between -1 and -2 MeV, here they range between -1.5 and -4.5 MeV. 
We attribute this difference mainly to the coupling with non collective phonons. As it was discussed alreday in Ref. 
[l3| . the energy shifts are mostly due to coupling with intermediate configurations including an octupole phonon. If 
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FIG. 8: (Color online) Feynman diagrams associated with particle [(a) and (b)] and hole [(c) and (d)] states. The shaded area 
denotes the RPA phonon and the wavy lines correspond to the interaction. See the text for further discussion. 
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FIG. 9: (Color online) Analysis of the (pi/2) ^ level density in ^''Ca produced by the coupling with 3 phonons. In the upper 
panels, all the calculated RPA 3^ phonons are included in the calculation. In the lower panels, only the lowest 3^ phonon 
state is instead taken into account. In the left panels, all possible intermediate single-particle configurations are included. In 
the right panel, the intermediate configurations are restricted to 8 3~. 
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FIG. 10: (Color online) The same as Fig. [9]in the case of the gg/2 level density. In the right panel, we show results obtained by 
taking into account only the coupling to specific intermediate configurations: P3/2 ^ 3~ (green curve), f5/2 ® 3~ (blue curve), 
and ir/2 ® 3~ (purple curve), respectively. 
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FIG. 11: (Color online) Energy of the seven neutron bound states lying close to the Fermi energy in ^''Ca, calculated in HF 
(left) and by taking into account the coupling to phonons (right). 
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we compare the theoretical results with experiment, we must probably conclude that a re-fitting of the effective force 
(SLy5 in the present case) is needed if this has to be used outside the mean-field framework. In fact, the HF-PVC 
results need a global upward shift in energy. 

1. Comparison with the experimental data in ^''Ca 




12 10 8 6 4 2 2 4 6 8 10 12 



Excitation energy in Ca [MeV] Excitation energy in Ca [MeV] 

FIG. 12: (Color online) The theoretical level densities pij (blue curves) are shown as a function of the excitation energy in 
^^Ca (hole states, left panels) and in ^^Ca (particle states, right panels). Except for the change of scale, the results are the 
same already shown in Fig. [T] They are compared with the experimental spectroscopic factors taken from [30l . [sij : these 
are represented in histogram form and also convoluted with Lorentzian functions of width 0.4 MeV (red curve). The vertical 
dotted line shows the one-neutron theoretical threshold energy obtained within HF-I-PVC (this energy is then the continuum 
threshold measured by setting at zero the energy of the renormalized lf7/2 state already displayed in Figs. [71 and [TT)| . 

The experimental single-particle strength of ''''Ca is obtained from ^°Ca(p, d)^^Ca pickup (for hole states) and 
''°Ca(d, p)^"'Ca stripping reactions (for particle states), by comparing the measured cross sections with Distorted 
Wave Born Approximation (DWBA) calculations performed with conventional assumptions. In particular, one usually 
assumes that the wavefunction of the transferred nucleon, (j)nij can be taken as an eigenfunction of a static mean field 
potential, by adjusting the depth of that potential so that the binding energy becomes equal to the experimental 
separation energy and the correct asymptotic dependence is guaranteed. The comparison with the level density 
obtained in a calculation like the present one, although not straightforward, is reasonable for levels which are well 
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FIG. 13: (Color online) The integrated level density as a function of excitation energy for ^^Ca (left panel) and *^Ca (right 
panel), respectively. The red curves correspond to the experimental data, while the blue curves are our HF+PVC calculation. 

described by the one-quasiparticle approximation. In the previous subsection, we have seen that this is indeed the 
case for several states close to the Fermi energy: for them, the single-peak associated virith a definite value of the 
number of nodes n, appearing in HF, persists. A diagonal, even perturbative, approximation for the mass operator 
is quite appropriate. However, for states characterized by a broad distribution in energy, when several values of n 
are mixed, the comparison with a simple DWBA calculation is likely to be less reliable (cf., e.g., the discussion in 
Ref. [3^). In principle, one should rather perform a direct theoretical calculation of the transfer cross section, using 
the wavefunctions that include many-body correlations. This goes beyond the scope of the current paper, and in the 
following we shall limit ourselves to a simple comparison with the spectroscopic factors reported in the experimental 
papers [s^, [Mj . Our results are comparable to those obtained in Ref. [33|] , where the distribution of single-particle 
strength in '^^Ca was calculated in a (discrete) quasiparticle-coupling model going beyond the diagonal approximation. 
The red histogram bars in Fig. 1121 show the experimentally determined spectroscopic factors, which are convoluted 
with Lorentzian functions having a width equal to 0.4 MeV to produce the red continuous lines. These can be 
compared with our theoretical level densities (blue continuous lines). The dotted vertical line shows the calculated 
threshold for one-neutron emission, which overestimates the experimental value by about 2 MeV. 

The total single-particle strength associated with the different quantum numbers l,j, obtained by integrating the 
level densities displayed in Fig. [12] up to E = 10 MeV is reported in Table IIIII One finds an overall satisfactory 
agreement between theory and experiment. The position of the centroid energies is reasonably well reproduced, 
except that in the case of the Si/2 strength, where the theoretical centroid energy is too low by about 2 MeV. In 
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Exp. 


Theory 




Exp. Theory 


d3/2 


0.88 


0.80 


17/2 


0.74 


0.66 


Sl/2 


0.84 


0.80 


Pl/2 


0.80 


0.81 


P3/2 


2.9 X 10" 


0.05 


P3/2 


0.73 


0.79 


d5/2 


0.73 


0.75 


d5/2 


0.11 


0.04 








l5/2 


0.88 


0.77 








§9/2 


0.28 


0.36 



TABLE III: Experimental spectroscopic factors Sij obtained from one-nucleon transfer reactions for hole and particle states in 
^®Ca and ^^Ca, compared to the integral of the theoretical level density performed up to an excitation energy of 10 MeV (cf. 
Fig. [131). 



general, theory tends still to underestimate the fragmentation of the single- particle strength: this occurs in particular 
for the strength (for particles and holes), and for the ic,/2 strength (for particles). This can be seen also from Fig. 
1131 where we show the cumulated experimental and theoretical strength distributions. The latter distributions tend 
to show a sharper increase. This can be attributed to several reasons. Among the possible ones, we point out that 
the present RPA calculation underestimates the experimental value of the collectivity of the low-lying 3j" phonon (cf. 
Table |ll]) and cannnot describe in detail the ISGQR strength distribution. If one could include its admixture with 
two particle- two hole configurations, these could shift part of the strength at lower energy and increase the effect of 
the coupling with the single-particle strength. A more fragmented ISGQR distribution would be in better agreement 
with experiment [13, H^; however, such a calculation would require to go beyond the formalism of the present work. 



B. Results for ^"**Pb 



Theory (RPA) Experiment 



Nucleus J" 


Energy 




Energy 






[MeV] 


[e^fm^-'] 


[MeV] 


[e^fm^-^] 


2"8pb 2+ 


5.12 


2.35 X 10^ 


4.09 


3.00 X 10^ 


3" 


3.49 


7.08 X 10^ 


2.62 


6.11 X lO'^ 


4+ 


5.69 


5.16 X 10^ 


4.32 


15.5 X lO'' 


5" 


4.49 


4.96 X 10* 


3.20 


4.47 X 10* 



TABLE IV: Energies and electromagnetic transition probabilities associated with the low-lying isoscalar collective states in 
^"^Pb. The experimental data are taken from Refs. ,28. ,361] . 

In Fig. [14] we provide an overall view of the calculated RPA multipole strength in ^"^Pb, while the energy and 
transition strength of the lowest states of each multipolarity are reported in Table IIVI The properties of the low- lying 
states are reproduced reasonably well by our calculation, with the partial exception of the transition probability 
associated with the 4+ state. In Fig. [151 we show the results of our systematic calculation of the level densities, 
displaying the outcome of the full HF-I-PVC calculation including 2+, 3~, 4+ and 5~ phonons (blue curve), as well as 
the results obtained by including only the 3~ phonons (red curve), in comparison with the HF results (black curve). 

The HF single-particle spectrum calculated for ^°^Pb is reported in Table |V] and illustrated in the left column of 
Fig. [16] and Fig. [TTl In the central column of Fig. [TBI and Fig. [ITjwe show the position of the main peaks obtained 
from the full HF-I-PVC calculation, whereas in the right column we include the experimental results. From an overall 
look at Fig. [151 we can notice that the quasiparticle picture (a single peak emerging from the PVC calculation, with 
shifted energy and renormalized integral with respect to HF) holds for most of the valence hole states, that is, for 
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FIG. 14: (Color online) The RPA strength functions in ^°*Pb obtained from our RPA calculations. The IS 2+ strength is 
reduced by a factor 5, while the IS 3~ strength is reduced by a factor 10. 



3pi/2, 2f5/2, 3p3/2, lii3/2 and lh9/2. A partial exception is constituted by the state 2f7/2, that acquires a double 
structure mainly due to the coupling with the 3j~ (ED ii3/2 configuration. The inclusion of PVC brings the relative 
position of the valence hole states in much better agreement with experiment. For particle states, the quasi-particle 
picture seems to be valid for 2gg/2, liii/2i 3d5/2, 4si/2, 2g7/2, 3d3/2 and, only to some extent, for lji5/2. The position 
of the 2gg/2, liii/2 and lji5/2 states is in good agreement with experiment, while the 3d5/2, lji5/2 and 2g7/2 orbitals 
lie too high in energy. 

The present calculation is similar to the one of Ref. [3| but the energy shifts are larger, as in the case of ^°Ca, 
due probably to the contribution of non-collective states. We can also make an overall comparison with the various 
results reported in Ref. [3] , by evaluating an average particle- hole gap defined as 



with 



{eh) 



Aoj = {cp) - (eh), 



Eu„occ(2j' + 1) ' 
Eocc(2j + l)enO- 
Eocc(2j + l) ' 



(22) 



(23) 



where the labels "unocc" and "occ" refer to the unoccupied and occupied valence shell, respectively. Starting from 
the HF value, 9.34 MeV, the HF-PVC value is reduced to 7.89 MeV and gets closer to the experimental value of 6.52 
MeV; the difference between the two values, namely -1.45 MeV, compares well with the values presented in Ref. Q, 
that range between -3.2 MeV and -1.1 MeV. 
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-9.1 
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nr 7 

2f2 


-12.1 




lg2 






lg| 


-27.5 




Ihf 


-12.8 




Ih^ 


-18.5 






-9.4 





TABLE V: Skyrme Hartree-Fock single-particle energies for Pb obtained with the force SLy5. Only states at negative energy 
are reported in the Table. 



The processes leading to the fragmentation of the single-particle strength for the orbitals lying far from the Fermi 
energy in ^"^^Pb have been already extensively discussed within the framework of more phenonienological studies 
however, for the convenience of the reader, in the following we present some details of the present calculation. 

In many cases, the strong broadening of the single-particle strength observed in Fig. [15] is caused mostly by the 
coupling with the ISGQR and ISGOR, due to the favourable matching with the difference between the relevant 
single-particle states. This is the case for the orbitals (lsi/2)~^ and (2si/2)~^) which couple to d3/2, d5/2 ® 2+ and to 
%|2^ h II ® 3~. The contribution of the 4+ strength is small, but not completely negligible. The main configurations 
contributing to the large broadening of the (3si/2)~^are f5/2, h I2 ^ 3~. In the case of (lpi/2)~^ and (2pi/2)~"'^, the 
2"*", 3~ and 4+ phonons give comparable contributions to the strength fragmentation, which is also quite large. In 
the case of (3pi/2) and (3p3/2)^^, there is no good match with the energy of available single-particle configurations, 
and this explains the small amount of fragmentation that characterizes these states. 

In the case of (ld3/2)~^, the 2+, 3^ and 4+ phonons give comparable contributions to the fragmentation. In 
the case of (2d3/2)^^, instead, the 3~ phonons play the most important role for the fragmentation: in fact, the 
relevant single-particle configurations are (3p3/2)~^, (2f5/2)~^j (2f7/2)~^ and (1119/2)"^ coupled with the low-lying 
3j~ state. A similar pattern holds for the spin-orbit partners, that is, in the case of (^5/2)^^, the 2+, 3~ and 4+ 
phonons give comparable contributions for the fragmentation but in the case of (2d5/2)^^, the 3^ phonons play the 
main role for the fragmentation (with some contribution arising from coupling with 4+ phonons). For (2d3/2)^^, the 
single-particle configurations involved are (3p3/2)~^, (215/2)""^, (2f7/2)^^ and (lhg/2)^^ (coupled with the 3j" state), 
as well as (lii3/2)^^- In the case of (115/2)""'^ and (lf7/2)~"'^, the 2+, 3^ and 4+ phonons give similar contributions 
to the strength fragmentation; however, the main configuration involved turns out to be (lii3/2)~^ ® 3~. As already 
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FIG. 15: (Color online) The same as Fig. [3 in the case of ^°*Pb. 
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mentioned, the state {2i-^/2)~^ is not affected much by the particle- vibration coupUng. In the case of the states (lg7/2) 
and (lhg/2)~^, the 3~ is the main responsible for the couplings; however, the fragmentation is rather small, and the 
energy shift is also small. In the case of (lg9/2)~"'^, the 2+, 3~ and 4+ phonons give comparable effects. In the case 
of (lhii/2)~^j the fragmentation is caused by the coupling with the configurations (3p3/2)^^, (2f5/2)~^i (2f7/2)~^ and 
(lhg/2)~^ <8> 4+. The state 2h2i/2 is a resonant state in the continuum: 2+ and 3~ give the main contributions to 
fragment its strength: lji5/2 (8> 2+, 3d5/2, 2g7/2: 2gg/2, liii/2 ® 3~ are the main states that produce the strength 
fragmentation. Also lji3/2 is a resonant state in the continuum: in this case, 3~ and 5~ are the most relevant phonons 
for the fragmentation of the strength: the main configurations are 2g7/2, 2g9/2, liii/2 €5 3~ and 3d3/2, 3d5/2, 2g7/2, 
2g9/2j liii/2 ^ 5~. Finally, in the case of the state lji5/2, the fragmentation is mainly caused by the coupling with the 
configurations liii/2, 2g9/2 ^ 3~. Once more, from considerations related to the matching of initial and intermediate 
energies, we expect that the low-lying 3^ state is the main contributor. 
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FIG. 16: (Color online) Neutron particle states in ^"^Pb. 
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FIG. 17: (Cilor online) Neutron hole states in ^"*Pb. 



The cumulative level densities of the various orbitals are compared in Fig. [18] to spectroscopic factors obtained from 
the ^°^Pb(^He, a) stripping reaction for hole states, and from ^"*Pb(d, p) reaction [33, for particle states. The 
integrated level density is compared to the sum of experimental spectroscopic factors in Table I VII The quasiparticle 
character of the orbitals lying close to the Fermi energy is a general result of our adopted theoretical framework. 
Our results are in fair overall agreement with the experimental findings from transfer reactions, which were able to 
locate most of the quasiparticle strength associated with the orbitals lying close to the Fermi energy. The quality 
of the agreement varies from one case to the other, and it is hard to decide whether this has to be attributed either 
to specific features of our model, or to deficiencies in the experimental extraction of the spectroscopic factors (as 
testified, e.g., by the fact that some of them exceed the maximum allowed value of one). Furthermore, we must 
recall that (e,e'p) experiments lead to much smaller spectroscopic factors, and that the relationship between the two 
kinds of experiments, as well as the reletive role of long- and short-range correlations, is a matter which continues 
to be actively debated [sol. l40l|. Last but not least, the very possibility of extracting a spectroscopic factor as a true 
observable, has been recently questioned [4l|, \^ . 
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FIG. 18: (Color online) The same as Fig. [T2]in the case of ^"^^Pb. The experimental data are taken from [stI. [3^. 
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TABLE VI: The same as Table [Till for 2°*Pb. 
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FIG. 19: (Color online) The same as Fig. [13] in the case of ^°*Pb. 



C. Results for 



In this subsection, we finally give results for ^''O, as an example of neutron-rich, weakly bound (S'„ = 4.1 MeV) 
nucleus. This nucleus is a doubly magic isotope, due to to the usual proton shell closure at Z=8 and to an "exotic" 
neutron shell gap appearing at N=16. The magicity of ■^'^0, had been suggested by theoretical studies ^43.-45] , and 
has been established by the measurement of the (unbound) ^^O ground-state and of its decay spectrum to ^^O, 
and by the extraction of the N=16 single-particle gap from the 23,24,25q ground-state energies [46l l47l |. As already 
mentioned, having a tool that allows studying weakly-bound systems by taking proper care of the continuum, is 
one of the main motivations of the present work. Interestingly enough, it has been suggested that nuclei in which 
the neutron separation energy becomes smaller than the proton separation energy are characterized by larger single- 
particle spectroscopic factors or, in other words, by more pure single-particle states. This is the feature emerging by 
the plot shown in Fig. 6 of Ref. and consequently, one of those that can be analyzed within our framework. We 
will come back to this point below. 

In Table IVlIl we provide the HF single-particle spectrum for neutrons. In Fig. [501 we illustrate our results for the 
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RPA strength functions. Experimental information, although scarce, is available. The main results are that (i) there 
should be no bound excited state [4§|, and (ii) the lowest excitation should be a 2+ state lying at 4.72 MeV In 
our RPA spectra, the lowest peak among those found for the chosen multipolarities is indeed a 2+ one, and its energy 
and electromagnetic transition probability are 3.4 MeV and 4.2 e^ fm'*. Experiment has also provided indications for 
the existence of a 1+ state lying at 5.33 MeV, but our calculations are limited to natural-parity states. In the case 
of this nucleus we compute the 1~ strength as well. To ensure that coupling with 1^ phonons docs not introduce 
any error associated with spurious strength associated with the translational mode, we follow the procedure already 
discussed in our previous paper [50] . We find a significant amount of dipole strength lying at energies somewhat below 
the usual (IS and IV) giant dipole resonances. 
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FIG. 20; (Color online) The RPA strength functions in obtained with the present continuum RPA. 



Nucleus hole states [MeV] 


particle states [MeV] 


-O Isi 


-39.2 


Idf -1.3 


2si 


-5.3 




Ipi 


-17.2 




Ipf 


-22.0 




Idf 


-7.5 





TABLE VII: Skyrme Hartree-Fock single-particle energies for O with the interaction SLy5. Hole and particle states at 
negative energies are displayed. 

In Fig. 1211 we display our results for the level density of ^'*0. Before entering into some detail, we discuss the main 
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emerging features and compare with what is known experimentaUy. In our calculation, the ld3/2 and 2si/2 states have 
a marked quasiparticle character, namely they are associated with a single narrow peak. It makes sense, therefore, 
to compare the experimental value of the gap with the HF and HF-PVC results for the energy difference between 
the ld3/2 and 2si/2 states, that is, 4.0 and 4.6 MeV respectively. The HF-PVC result is in good agreement with the 
experimental value of 4.86 MeV. In heavy nuclei, as a rule, the PVC shrinks the single-particle gap and increases the 
effective mass (cf. the previous subsection) , but this is not the case in light nuclei due to the specific effect of having 
only low angular momentum occupied states (as already noticed in Ref. [l3|). In the present case, while the PVC 
pushes the ld3/2 orbital closer to the Fermi energy (-2.5 MeV compared to the HF value of -1.3 MeV), the 2si/2 hole 
state is pushed further from it (-7.1 MeV compared to the HF value of -5.3 MeV). 

The peak energies of the other orbitals obtained by using HF-PVC (HF) read 2.4 MeV (4.3 MeV) for lf7/2, and 
-8.3 MeV (-7.5 MeV) for ld5/2. The net effect of PVC is a shift down of the states. The absolute value of the energies 
is expected to depend on the choice of the effective force. Skyrme forces, as other mean-field frameworks, tend to 
predict larger binding in light neutron-rich nuclei as compared with the experimental findings, as it is clearly testified 
by the fact that ^^O turns out to be bound in many of these models. In the present case, the ld3/2 state is bound 
while it should be a resonant state. We can nonetheless look at relative energy differences. The known states in ^'^O 
taken from Ref. [U are, in addition to the 1/2+ ground-state, a 5/2+ state at 2.79 MeV and a 3/2+ at 4.04 McV 
(leaving aside the state at 5.34 whose character is not clear, being either 3/2~ or 7/2"). These are states that can 
decay to the ^^O ground state. In our calculation we can identify states below the energy threshold for this kind of 
decay: in particular the first 5/2+ state lies at 1.2 MeV in our calculation. 




FIG. 21: (Color online) The same as Fig. [7] in the case of ^*0. In this case not only 2+, 3 and 4+ phonons but 1 phonons 
as well are considered in the case of the full PVC result (blue curve). 

We now discuss, the couplings that produce fragmentation of most of the single-particle strength distributions. 
At variance with the state 2si/2, the state lsi/2 is strongly fragmented. This fragmentation is chiefly caused by the 
configurations IdT/j 2+ and (lpi/2)~"^, (lP3/2)~"'^ By inspecting the energy difference, we can assume that the 
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24q 

Holes Particles 

r Exp. Theory r Exp. Theory 
Si/2 — 0.81 d3/2 — 0.83 
ds/a — 0.78 

TABLE VIII: The same as Table HI] for ^^O. 

IVGQR and IVGDR play the main role for the fragmentation. The state (lpi/2)^^ is fragmented due to the coupling 
with the configuration (ld5/2)^^ ® 3^: we expect that the ISGOR and IVGOR play the main role by considering 
the energy matching. (2si/2)~^ ® 1~ can also contribute to the fragmentation, yet to a minor extent. The main 
configuration giving rise to the fragmentation of the state (lp3/2)~^ is the configuration (ld5/2)~^ ® 1~. Here, due to 
the energy difference between the hole states (lp3/2)~^ and (ld5/2)~"^, the dipole excitations around 14 MeV do play 
the main role. In the case of the fragmentation of the (ld5/2)~^ state, the configuration (2si/2)~^ ® 2+ is the most 
important one, and the low-lying 2+ state at 3.4 MeV is the most relevant. Finally, the energy shift of the state lf7/2 
is mostly caused by the coupling with the configuration ld3/2 ® 3~. 

Last but not least, we come back to the point raised above, namely that the (quasi-particle-like) states around the 
Fermi energy, ld3/2 and 2si/2, are quite pure (see also Table |VlTT|. In our calculation, there are not optimal energy 
matching of those states with other configurations due to the scarcity of low-lying collective excitations and large 
gaps between single-particle state. More generally, in our calculations the coupling of neutron states is mainly with 
the proton component of the phonon states (due to the dominance of the neutron-proton interaction). Therefore, 
neutron states on neutron-rich nuclei are expected to be more pure because proton excitations are pushed at higher 
energy as the neutron excess increases. Recently, other calculations of the spectroscopic factors in light nuclei within 
the coupled-cluster approach have become available in Ref. (52i] (see also the critical discussion in Ref. [ssj). 

IV. SUMMARY 

The idea that single-particle strength is not systematically pure in atomic nuclei, and that coupling with other 
degrees of freedom is quite relevant, is an old idea in nuclear physics. Phenomenological calculations based on 
particle-vibration coupling (PVC) for spherical nuclei have been performed for several deacades, and they have been 
quite instrumental to point out some of the limitations of the pure mean- field approach (like, e.g., the enhancement of 
the effective mass around the Fermi energy) . Microscopic PVC calculations based on the consistent use of an effective 
Hamiltonian, have become available only recently. 

None of the mentioned calculations, to our knowledge, takes proper care of the continuum. In our work, for the 
first time, we have implemented a scheme based on coordinate-space representation in which the single-particle states, 
the vibrations, and their coupling are calculated with proper inclusion of the continuum. This is of special interest if 
weakly-bound nuclei close to the drip lines are to be studied. However, also in well-bound nuclei the present approach 
present advantages in the sense that resonant states can be properly studied. Transfer to the continuum has been the 
subject of several experimental studies. 

In stable nuclei we obtain results that are in overall agreement with previous studies. We can, at the same time, 
better describe the fragmentation of the single-particle strength. The shifts of the single-particle states around the 
Fermi energy, with respect to the HF values, are relatively large (in keeping also with the fact that we cannot restrict 
our coupling to collective states only). We obtain an overall agreement with experiment in ^°®Pb, while in the case of 
^''Ca our results point to the need of re-fitting the Skyrme interaction that has been devised to work at the mean- field 
level and not beyond it. 

We have also applied our model to a neutron-rich nucleus, namely '^^O. This is a double-magic nucleus, and there 
are few low-lying states. Because of this, and also since the neutron states energy would be more affected by coupling 
with protons, the neutron single-particle strength around the Fermi energy is quite pure (i.e., spectroscopic factors are 
rather close to one). While this is in agreement with some experimental findings, certainly more detailed spectroscopic 
studies are needed to extract a global trend and a firm understanding. 
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Appendix A: Hartree-Fock Green's function 

As it is stated in the main text, it is necessary to use the causal Green's function for the Dyson equation because this 
equation is based on the use of the Wick's theorem, and the Wick's theorem applies only to time-ordered products. 
On the other hand, the continuum HF Green's function is given in the form of a retarded function. So we need to 
compute the causal Green's function starting from the retarded Green's function in order to use the continuum HF 
Green's function in the Dyson equation. 

The causal HF Green's function is defined by 

iG^'{rat,r'a't') = (*o|T{^(rai)?/;t(r'crY)}|*o) (Al) 
- B(t - i')(*olV'(rffi)^^(r-Vt')l*o) - ^(i' - t){^^\i,\r' a't')i,{rut)\^^). 

(A2) 

The retarded Green's function is instead defined by 

iGl^{r<Tt,r'a't') = 9{t-t'){'fo\{'4^{rat),'ti;^{r'a't')}\'io). (A3) 
From these two definitions, and using 9{t' — t) = 1 — 9{t — t'), we can find that 

iG^ircrt, r'a't') = iG^{rcrt, r'a't') ~ ^ e-^^'^^''*">(l)l{r'a')Mrcr) (A4) 



The Fourier transform of Eq. (|A4p is expressed as 



G^ircT,r'a'-u;) = (ra, rV; c.) + J] ^ ^J^T^^.^ramr'a') (A5) 

(in the limit 77 ^ 0) (A6) 
^ G^ {ra,r' a' ; uj) + 2mJ2Si^ ' eh)Mr<j)cbUr''j') (A7) 

h 

The continuum HF Green's function Gq.ij given by Eq. ([T5|) is regular for the complex energy E. So the retarded 
Green's function with a smearing width 77 is expressed as 

G^y irr' , w) = Go,ij {rr' ; w zt?) , (AS) 

where uj is the real part of the complex energy. Then the continuum causal Green's function can be expressed by 

Go,ij{rr',uj)^Go,ij{rr';uj + iT])+ ^ — — ^^^^j— ^0n;,i;,j„ (r)0„^/^j„ (r'). (A9) 



Appendix B: Unperturbed response function 

In general, the RPA theory can be formulated in two ways. One is based on the causal function, while another 
one on the retarded function. Both formulations give the same results for the physical part of the spectrum, namely 
for positive excitation energy. However, a complete RPA basis must include negative energy states and the two 
aforementioned formulations are different for negative energies. 

In order to construct the self-energy function of Eq. ([7]), we need to consider the RPA response function not only 
in the positive energy domain but also in the negative energy domain due to the required energy integration. So 
we cannot use the retarded RPA response function for the self-energy function. The RPA equation for the response 
function is given by i? = Ro + RqkR in any of the representations, Rq being the unperturbed response function, and R 



27 



the RPA response function. In order to obtain the causal RPA response function by solving this equation, the causal 
unperturbed response function should be used. 

Normally the continuum (Q)RPA is formulated by using retarded functions within the linear response theory. The 
continuum HF(B) Green's function is used to build the (retarded) unperturbed response function [13, HI]- I* 
therefore necessary to know how to convert the retarded unperturbed response function to the causal function in the 
continuum RPA formalism. We show it in the present Appendix. 

The causal and retarded response function in RPA are defined by 

iR^irt,r't') = {^o\nSpirt)Sp{r't'mo) (Bl) 
= e{t - t'){^o\Sp{rt)Sp{r't')\-^o} + e{t' - t){^^\5p{r't')5p{rt)\mo) (B2) 
iR^irt,r't') = eit ~ t'){^o\[6pirt),Spir't'mo). (B3) 

respectively (here j^I'o) is the HP ground state). From these definitions, one can find the relation between them as 
follows, 

iR^{rt,r't') = iR^{rt,r't') + {^>o\5p{rt)5p{r't')\^o) (B4) 

hh' 

The Fourier transformation of the latter equation gives 

R^irr'-u) = R^irr'-^)~Y.l ^ <l>l{r)4>h'{r)<t>l,{r')Mr> (B6) 

(w - eh' + ehY + ?7 

(in the limit of ry 0) (B7) 

^ <(rr';c^)-2^*^ (5(0. -ev+ 6,00;; (r),^„.(r)</.*,(r').^h(r'), (B8) 

hh' 

where RQ{rr';uj) can be expressed by means of the retarded HP Green's function as 

<(rr';L.) = ^ 0;;(r)G^(rr'; + e;,)(/.;,(r') + (/.;;(r')Gf(rr'; -c. + e,)0,(r). (B9) 

h 

In the continuum RPA formalism, the continuum HP Green's function is used as Gq in Eq. (jB9l) . 

Appendix C: Spectral representation of the Green's function and the response function 

Here we show the spectral representations of the HP Green's function and the RPA response function (both causal 
and retarded). The difference will appear in the sign of the imaginary part 77. Actually this is very important to 
obtain the proper self-energy function by using the contour integration, because this sign differencee produce changes 
in the position of the poles of the Green's function and of the response function on the complex energy plane (this fact 
is connected with the fact that the Wick's theorem can be applied only for the causal function, as already mentioned). 

^-^ UJ ~ eh ~ iri ^-^ id — eri + in 

h P f I 

G-ir.,rW;u) Mr^mrW) ^y^Mr^^jrW) ^^^^ 

^ LU - eh +iri ^ LU - e„ + iri 

h P f I 



Rc{rr';uj) =^ 



(0|p(r)|^)(Hp(O|0) (0|/3(rO|^)(H/5(r-)|0) 
UJ — + iri UJ + Eu — irj 

(G3) 

R irr'-uj) V (0|p(r-)|^)(Hp(rO|0) (0|p(rO|^)(H/5(r)|0) 
^■^ ' > 2^ i^_E„+ir] uj + E^ + iT] ^ ' 
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Appendix D: Spectral representation of the self-energy function 

The self-energy function in the space-time representation is defined by Eq. (|4]). If we insert the HF and the RPA 
resuhs in this definition, then the self-energy function can be expressed as 

S](riaiii,r2(T2t2) = K{ri)Go{ri(Titi,r2(J2t2)K{r2)iR{ritir2t2) (Dl) 
= ]e{h - t2) e-*(^'^+^''^^*^~*^)K(ri)<5p.(ri)0p(riai)^;(r2a2)(5p:(r2)At(r2) 

--Oit2 - h) V e+'(^-^'')(*^-*^)A^(ri)<5p:(ri)</.;;(riai),^,(r2a2)5p.(r2)«:(r2) (D2) 
t ^ — ' 

The Fouriter transform of Eq. (jD2p gives 

M r'y)Sptir)Kir)^lir'a')Sp,ir')Kir') 

Ul - Ch + Ey - IT] 

h,v 

^ ^ c^j,[ra)5p.[r)n{r)(^;[r'a')5pl{r')n{r') 



Appendix E: Residual interaction within the Landau-Migdal approximation 



Here we show the explicit expression of the residual interaction within the so-called Landau-Migdal approximation. 
This residual force is used in the self-energy function [Eq. (fT3|) ]. 



Kgg'(r-) 



dh„ , , S^E 
-(r) = — — — 

dpqdpg' 

to. 



dpq' 



where q,p 



(9 -9') f(l-2;o) 



-12^' 



(7 + 2)(7+l)(l + f )-(.X3 + ^)(2 + 47^+7(7-l)E( V 



Pq 



+ -{tl{l-Xi)+it2{l+X2))kl 



= < 



{q^q') io(l + y) 



-12^' 



(7 + 2)(7+l)(l + f )-(x3 + ^)(27 + 7(7-l)E(^ 
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